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We prove the existence of spatially localised and time-periodic solutions in general nonlinear Hamil¬ 
tonian Klein-Gordon lattice systems. Like normal modes, these localised solutions are characterised 
by collective oscillations at the lattice sites with a uniform time-dependence, an approximation that 
is valid for lattice states near the end points of the linear spectral band. The proof of existence 
uses the comparison principle for differential equations to demonstrate that at each lattice site every 
half of the fundamental period of oscillations, contributing to localised solutions of the nonlinear 
lattice, is sandwiched between two oscillatory states of auxiliary linear equations. For soft (hard) 
on-site potentials the allowed frequencies of the in-phase (out-of-phase) localised periodic solutions 
lie below (above) the lower (upper) value of the linear spectrum of phonon frequencies. A numerical 
analysis of the solution behaviour verifies that localised solutions with frequencies close to the edges 
of the the linear phonon spectrum exhibit (virtually) uniform time dependence. 

PACS numbers: 05.45.-a, 63.20.Pw, 45.05.-l-x, 63.20.Ry 


I. INTRODUCTION 

Ever since their occurrence was proposed in [1], intrinsic localised modes (ILMs) or discrete breathers in nonlinear 
lattices have attracted continued interest, not least due to the important role they play in many physical realms where 
features of localisation in systems of coupled oscillators are involved (for a review see [5] and references therein), [3]-[^. 
Proofs of existence and nonexistence of breathers, as spatially localised and time-periodically varying solutions, were 
provided in [6]-[14] using different techniques. The exponential stability of breathers was proven in [15]. Analytical 
and numerical methods have been developed to continue breather solutions in conservative and dissipative systems 
starting from the anti-integrable limit m-m- 

During recent years the existence of ILMs (breathers) has been verihed in a number of experiments in various con¬ 
texts including micro-mechanical cantilever arrays [22] . arrays of coupled Josephson junctions [23], antiferromagnetic 
chains [33], coupled optical wave guides [l^, Bose-Einstein condensates in optical lattices [35], in coupled torsion 
pendula [37], electrical transmission lines [35], and granular crystals [33]. Regarding their creation mechanism in 
conservative systems, modulational instability (MI) provides the route to the formation of breathers originating from 
an initially spatially homogeneous state imposed to (weak) perturbations. To be precise, the MI of band edge plane 
waves triggers an inherent instability leading to the formation of a spatially localised state [30] . 

Lattice discreteness and nonlinearity are the prerequisites for the ability of breathers to store energy without 
dispersing it to the environs. Discreteness serves for the existence of an upper boundary to the linear phonon spectrum 
making it possible that not only the fundamental breather frequency but also all of its higher harmonics lie outside 
of the phonon spectrum. Nonlinearity allows to tune the frequency of (anharmonic) oscillators in dependence of their 
amplitude (or action) such that resonances between harmonics of the breather frequency, ujb, and the frequencies, 
LOph, of linear oscillations, i.e. phonons, can be avoided. That is, the non-resonance conditions 

mwf, 7 ^ OJph, Vm G Z (1) 


need to be satisfied. 

In this work we prove the existence of spatially localised time-periodic solutions for which the small amplitude 
oscillations at the lattice sites are characterised by a uniform time-dependent functional form. The paper is organised 
as follows: In the next section we introduce general nonlinear Klein-Gordon lattice system and discuss the features 
of time-reversibility. Subsequently the theorems together with their proofs of the existence of time-periodic and 
spatially localised solutions in KG systems with soft and hard on-site potentials are presented. It follows a section 
on the numerical analysis of the solution behaviour of lattice states with frequencies close to the edges of the linear 
phonon band. Finally we summarise our results. 
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II. NONLINEAR KLEIN-GORDON LATTICES 

We study the dynamics of general one-dimensional infinite nonlinear Klein-Gordon (KG) lattice systems given by 
the following system of equations 


(Pqn 

dt'^ 


U (^Qn) “f Qn—1)5 Tl ^ ^ 


( 2 ) 


where the prime ' stands for the derivative with respect to the latter being the coordinate of the oscillator at site n 
evolving in an anharmonic on-site potential U{qn)- Each oscillator is coupled to its next neighbours and the strength 
of the interaction is regulated by the value of the parameter n. This system has a Hamiltonian structure related to 
the energy 



(3) 


and it is time-reversible with respect to the involution p i— >■ —p. 

The on-site potential U is analytic and is assumed to have the following properties: 

17(0) = C7'(0) =0, 17"(0)>0. (4) 

Furthermore, we assume that U{q) = U{—q). In what follows we differentiate between soft on-site potentials and 
hard on-site potentials. For the former (latter) the oscillation frequency of an oscillator decreases (increases) with 
increasing oscillation amplitude. A soft potential possesses (at least) two inflection points at —qi < 0 and Qi > 0 and 
it holds that 

U\-q^ < <? < 0) < 0, C/'(0 <q<qi)>G (5) 

U"{±qi) = 0, U"{-q, <q<qi)>0. (6) 


We remark that U{q) can have more than two inflection points (an example is the periodic potential U{q) = — cos{q)). 
In the frame of the current study we are only interested in motion between the inflection points adjacent to the 
minimum of U{q) at g = 0. 

Hard on-site potentials are, in addition to the assumptions in Q, characterised in their entire range of definition 
by 


U'{q < 0) < 0, U'{q > 0) > 0, U''{q) > 0. 
The linearised system admits plane wave (phonon) solutions with frequencies 

Wp?j(fc) = wg+4Ksin^ , fce[-7r,7r] 


and ClIq = 17" (0). 

We study solutions to the system (§ that are time-periodic satisfying 

qn{t + Tb) = Qnit) 


(7) 

( 8 ) 


(9) 


with period Tf, = 21^jLOb- Moreover, these periodic solutions are supposed to be localised about a single site, say n = 0, 
exhibiting the reversion symmetry qnit') = q_„(t). 

In compact notation the canonical equations can be expressed as 

X = F(X), (10) 

where X = (Ar„)„gz is the infinite sequence of pairs A„ = (p„, of canonically conjugate momentum and coordinate 
variables and the vector field F = (F„)„gz is determined by 


Fn = {-dHjdqn,dHldpn). 


( 11 ) 


In particular, due to the finite radius of the interaction between the oscillators (nearest-neighbour coupling) the r.h.s. 
F of the system (101 is of the form Fn = F„(A„_i, A„, A„+i). 
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Introducing the involution 



the time-reversibility of the system (101 is reflected in 


pF{X) 


-F{p{X)). 


( 12 ) 


(13) 


The forthcoming results in relation to the existence of localised periodic solutions to 
following statement. 


(i) are supported by the 


Lemma: Let X(t) = (X„(f))„gz be a solution to the time-reversible system with the reversible vector field (ll\ l 
and with initial conditions X(0) = (X„(0))„gz such that Pn(0) 7 ^ 0 and q-niO) = 0 for n G Z and there is some T > 0 
such that pX(0) = X(T). Then on the interval [0,T] one has pX(t) = X{T — t), i.e. X(t) is reversible on [0,r]. 


Proof: The solution to X = F{X) with initial conditions X(0) for t G [0,r] is expressed as 

X(t)=X(0)-f f dTF{X{T)). 

Jo 

Then, using the reversibility properties pX{0) = X{T) and pF{X) = —F{pX), one gets 

pXit) = pXiO) + p [ drFiXir)) 

Jo 

= XiT) + [ drpFiXir)) 

Jo 

= X{T) - f drF{pX{T)) 

Jo 

= X{T-t) 


(14) 


and the proof is complete. 

The next section deals with the existence of localised solutions. 


□ 


III. EXISTENCE OF LOCALISED SOLUTIONS 

In the following we prove the existence of spatially localised and time-periodic solutions for the system where 
the oscillators perform collective motion of the form qn{t) = 4>nf{t) with the sequence {(j)n)nGi determining the 
amplitudes of periodic oscillations, the latter being represented by a periodic function f{t) = f{t T^) of period 
T = 2TTfuJii [3T]. Localisation requires lim|„|_,,_|_oo Wfin — 4>n-i\\i°° = 0- These solutions, like normal modes, are 
understood as a vibration in unison of the system, i.e. all units of the system perform synchronous oscillations so that 
all the oscillators pass through their extreme values simultaneously [32],[33]. The latter feature is also common to 
discrete breather patterns arising in nonlinear lattice systems, where the oscillators perform in-phase and/or out-of- 
phase single-frequency periodic motion with respect to a reference oscillator [18j.|35j. By comparison to the uniform 
time-dependence of NNMs, in general for discrete breathers the functional form of the time-dependence changes from 
site to site. Nevertheless, for small amplitude localised solutions in infinite lattices with frequencies of oscillations, 
Wf,, near the end points of the spectral phonon band (i.e. close to the lowest (highest) normal mode frequency in 
systems with soft (hard) on-site potentials, i.e. Wh ^ wq {tOb ^ -^/wq -I- 4k)) solutions are well approximated in the 
form g„(t) = finfit) [S]. In [S],[II], using a center manifold technique, it is proven that small amplitude breather 
solutions in one-dimensional infinite nonlinear lattice systems, possessing a band formed by the frequencies of the 
linear spectrum, exist and are of the form qn(t) = 6 „/(f) where f{t) is a periodic harmonic function in time and 
the decaying amplitudes bn constitute an exponentially localised shape as |n| —>■ 00 . Such discrete breathers of small 
amplitude emerge from a bifurcation of a band edge plane wave. In the limit of their amplitude approaching zero the 
degree of localisation (localisation length) goes to zero. Furthermore, approximations of small amplitude breathers in 
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KG lattice systems near the end points of the spectral band by solitons of discrete nonlinear Schrodinger equations 
can be obtained by means of a resonant normal form theorem [84] . We represent here a proof of the existence of 
small amplitude time-periodic and spatially localised solutions in infinite one-dimensional KG lattices by means of an 
alternative technique. 

Based on the work in mm. we make the Assumption that lattice states with frequencies close to the phonon 
band edges are well approximated by 9n(t) = (f’nfit) with n S Z, and state the main results of the paper showing 
that the oscillators perform in-phase (out-of-phase) periodic motions in systems with a soft (hard) on-site potential 
with decaying amplitudes (/)|„| > ±n S N forming a localised pattern. 

Theorem 1: Consider system ^ with a soft on-site potential. Suppose the solutions corresponding to lattice states 
with frequencies close to the lower edge of the linear phonon hand are of the form 

qn{t) = (fnfit), n e Z, G M, (15) 

fulfilling |9n(t)| £, < Qi fort G K where ^ and qt denote the maximal amplitude of motion and the (positive) inflection 

point of the soft potential respectively. Provided the inequality 

U"{0 > nA, (16) 

where 

A f (^n+1 ‘^4'n 4” 4*n—l) (0n-t-2 4” ^n) 

A = max < --- --- 

neZ+ ( (pn (pri+l 

is satisfied then there exist localised periodic solutions {qn{tTt), qn{tTb)) = (<?«(!), 9n(t)), so that the oscillators 
perform in-phase motion, i.e. sign(g„) = sign(g„+i), given by with 0 < 4>n+i < 4>n for n > 0 and (pn > 4>n-i > 0 for 
n < 0. f is a periodic smooth function of time t with period Tf, = 27r/wh and frequency ujb satisfying 



ujb < Wo- 


(18) 


For time-reversible solutions f{t) is determined by f(t) = cos(uJbt). 


Proof: W.l.o.g. the initial conditions are chosen such that gn(0) = 0, and 4^(0) 0, say qn{0) > 0, Vn, and 

g„+i(0) > gn(0) > 0 for n < 0, and qn{0) > 9n-i-i(0) > 0 for n > 0. Then due to continuity there exists some t, > 0 
so that during the interval (0, G] the following order relations are satisfied 

0 < qn+i{t) < qnit) , n>0, (19) 

0 < qn{t) < q„+iit ), n < 0. (20) 

We define the difference variable between the coordinates and velocities respectively at sites n and n 4- 1 as follows 


Ag„(t) = a[qn{t) - g„+i(t)], Aqn{t) = a[qn{t) - g„+i(t)], 


( 21 ) 


where a = 1 {a = —1) for n > 0 (n < 0). Thus by definition Aqn{t) > 0 on 

Due to reversion symmetry q^ = q-n it suffices to consider n > 0, viz. a = 1. The time evolution of the difference 
variables Aqn{t) is determined by the following equation 


d?Aq„ 

df^ 


— — \U'{qn) — U'{qnj^i)] 4- K(Aq„+i — 2Ag„ -|- Aq^-i). 


Utilising that for qu > qi and |g„|„gz < ^ < it holds that 

U"mqk - qi) > U'{qk) - U'(qi) > U"{0{qk - qi) > 0 
enables us to bound the evolution of the difference variable Aq on from above as 

dfAqn 




< -U'fCjAqn 4- k(A9„+i - 2Aq„ -|- Ag„_i) , 


( 22 ) 


(23) 


(24) 







and from below as 
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d‘^Aq„ 

(W^ 


— ^qAqji 2Aqyi + Aqji_i'). 


Substituting g„(t) = 'pnf{t) into Eqs. (24) and (25) gives 


<Pl 

dip' 


<-[u'\0- 


{.4>n+l — ‘P‘4’n + (t>n-l) ~ (</’n+2 ~ ‘P‘4’n+1 + 4‘ri) 


and 


fl 

dt'^ 


> — \ UJq — K 


^n+1 


(0n+l ‘P‘4^n “b 4^n—l) (0n+2 20^_|_i + (ffn) 

4’n 4^n+l 


f <-{U"{0-^A)f, 


^ / > - (wo - kB) f. 


with A given in 0 and B is determined by 

(0n+l ‘P‘4^n “b 4^n—l) (0n+2 ‘^^n+l “b ' 


B = min 

nGZ+ 


i^n+1 


(25) 


(26) 


(27) 


(28) 


In compliance with the order relations (19) the amplitudes <()„ form a localised sequence 0 < (j)n+i < <Pn for n > 0. 

Therefore, by the comparison principle for differential equations, for given initial conditions, f(t) and df{t)/dt = f{t) 
are bounded from above and below by the solution of 


and 


respectively, provided G{t) > 0 and g{t) > 0 and we introduced 


d'^G 

(29) 


d^g .2 

(30) 

1 we introduced 



(31) 

(Ilf = UJq — kB. 

(32) 


By the hypothesis (16) it is assured that fig G K. Furthermore, since wq > fb"(C) it follows that (Ds S ffi and so 
there exist oscillatory solutions t o and (301. 

The solution to Eq. (29) and (pOlwith initial conditions (/(O) = 0,/(O) > 0) (complying with <7^(0) = 0 and 
q„( 0 ) 7 b 0 for all n) is given by 


and 


G{t) 

/(O) . . 

(33) 

Git) 

= /(O) cos(Dgt) 

(34) 

git) 

/(O) ■ ^ 

= - 

OJg 

(35) 

git) 

= /(O)cos(wgt) 

(36) 


respectively, where G{t) > 0 for 0 < t < Tr/Clg and g{t) > 0 for 0 < t < Tr/wg. 

Notice that dPf{t)/dt^ < 0 on (0,t*), that is, the acceleration stays negative. Due to the relations / > 0 on (0,<») 
in conjunction with the lower bound g{t) < f{t) the order relations as given in (19),(20) are at least maintained on 
the interval ( 0 , 7 r/wg). Moreover, f{t) is bound to grow monotonically at least during the interval (0, 7 r/( 2 a;s)] and 
attains a least maximal value f{Qi)/Cjs- Furthermore, f{t) cannot return to zero before t = tt/ws- 

On the other hand, from the upper bound f(t) < G{t) one infers that f{t) can attain an absolute maximal value 
/( 0 )/Os but not before t = tt/{2(1,s) and f{t) is bound to return to zero not later than t = 7 r/i 2 g. 

Therefore, there must be a S [tt/ws, tt/IIs] such that /(t*) = /(O) = 0 and /(t*) = —/(O). 
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Hence one has qn{t*) = 9n(0) and qn{t*) = —(7ri(0), and by the Lemma above their time evolution is given in terms 
of reversible solutions obeying 


qn{t*/2 + t) = qniU/2-t), (37) 

qniU/‘2 + t) = -g„(t*/2 - t)), (38) 

with 0 < t < t*/2 and at the turning point, t = t*/2, the coordinates qn assume simultaneously their respective 
maxima while the ()„ are all zero. Conclusively, on the interval [0, t*] the oscillators evolve through half a cycle of 
periodic in-phase motion, i.e. sign(g„) = sign(q„+i) and sign(g„) = sign(()„+i). At t = t* one has g„(t*) = qn(0) = 0 
and quit) = —qn{0) so that by symmetry during the interval [t*,2t*] the oscillators repeat their behaviour exhibited 
during the initial interval [0, <»] but this time with reversed sign of their coordinates and velocities and hence, at 
t = 2t* one period of localised periodic oscillations is completed. Repeating the process results in time-periodic and 
spatially localised motions of the oscillators. 

By taking the values of the coordinates at the (first) turning point t = t*/2 together with ()„ = 0 as (new) initial 
conditions one generates time-reversible solutions possessing the symmetry quit) = qn{—t) and —qnif) = 9n(—t) so 
that the order relations 


\qn+i(t)\ < \qnit)\, n>0, 
\qn+i{t)\ > \qn{t)\, n<0 


(39) 

(40) 


are true for t gM.. As lim|„|_^oo q-n = ^ the tails of the breather solutions behave linearly so that their time-dependence 
is characterised by a single-frequency harmonic function. The latter for time-reversible solutions is determined by 
f{t) = cos{uJbt). 

Furthermore, the amplitudes in the exponentially localised tails of the breathers can be represented as 4>\n\ = 
with 0 < A < 1, ztn G N and AT > 0 S K and when substituted in (17) one gets 


-)- (^n-l-2 ^^n+1 “t” ^n) 

4^n+l 


= K{X+^-2 


> 0 for 0 < A < 1. 


Then concerning the frequencies of oscillations one finds from (32) that 


UJb < 


{4^n+l 2i(l)>^ -\~ —l) (^An-t-2 2^^n+l “t” 

4^n 0n-t-l 


< Wo 


(41) 


(42) 


and thus the breather frequenci es li e belo w th e lower edge of the linear band. 

Note that the inequalities in (39) and (40) are strict except at those moments in time, tfcgz, when the oscillators 
pass simultaneously through zero coordinate, corresponding to the minimum position of the on-site potential, i.e. 
qn{tk) = 0 and the proof is complete. 

□ 


Remark: The relations 


{i’n+l — ‘2‘4>n + (f’n-l) ~ {4'n+2 ~ ‘2'4>n+l + 4>n) 


I t/"(f) — K max 
nez,+ 


— <pn+l 


< l^b 


^ 2 • f “f 4^n—l) (0n-t-2 -)- ((>n) 

< \ LO^ — K mm < -^^- 

y ^ ™ez+ ( 


f>n — (Pn+1 


(43) 


entail that 


lim lim uit = ojq , C G M. (44) 

This is in compliance with the fact that breathers of (arbitrarily) small amplitude exists in one-dimensional infinite 
anharmonic lattices near the (lower) edge of the linear band of phonons. In fact, small amplitude breathers emerge 
from band edge plane waves by means of a tangent bifurcation [lOj . Conversely, in the limit —>■ (7 with C G IR, i.e. 
when the localisation strength approaches zero, the breather is transformed into a spatially extended solution, viz. a 
plane wave. 
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Moreover, the condition (43) involving the difference between the (local) curvatures ((/>n+i — 2(/)„ + 4>n-i) and 
{ 4 'n +2 — ‘2‘4‘n+i + 4>n) imposes, in dependence of the maximal amplitude of motion,^, and the coupling strength k, a 
constraint on the (local) degree of localisation. 


The next Theorem establishes the existence of out-of-phase localised solutions for motions in hard on-site potentials. 

Theorem 2: Consider system ^ with a hard on-site potential. Suppose the solutions corresponding to lattice states 
with frequencies close to the upper edge of the linear phonon hand are of the form 

qn(t) = {-i)^(j)nf(t), n e Z, (()„ e K, (45) 

fulfilling \qn(t)\ < ry fort € M where rj denotes the maximal amplitude of motion. There exist localised periodic solutions 
iqn{t+Tb),qn{t-\-Ti,)) = iqn(t),qn{t)), so that the oscillators perform out-of-phase motion, i.e. sign(g„) = — sign(q„+i) 
given by ( |'^5P and 0 < 4>n+i < 4>n for n > 0 and 4>n > 'Pn-i > 0 for n < 0. f is a periodic smooth function of time t 
with period Tf, = 27: jiUb and frequency ojb satisfying 

ojb > -I- 4k. (46) 


Proof: W.l.o.g. the initial conditions satisfy 


g„(0) = 0, 0< (-l)"+i(?„+i(0) < (-l)"g„(0), n > 0, (47) 

g„(0) = 0, 0 < (-l)”g„(0) < (-l)”+ig„+i(0), n < 0. (48) 

Then due to continuity there exists some > 0 so that during the interval (0, t*] the oscillators perform out-of phase 
motion and the following order relations are satisfied 

0 < (-l)”+^g„+i(t) < (-l)”(y„(t), n > 0, (49) 

0 < (-1)”9„(<) < (-l)"+ig„+i(t), n<0. (50) 


We proceed as in the previous case for soft on-site potentials by introducing the difference variable between coordi¬ 
nates as in (211. The time evolution of the difference variable is determined by an equation similar to (221 and using 
that for qk > qi it holds that 


U^'iOiqk - qi) > Ufqk) - U\qi) > t/"(0)(gfc - g,) > 0 (51) 

we derive on [0,t*] bounds for evolution of the difference variable Aq^ft) an upper bound and lower bound as 

dfAq„ 


and 




dfAqn 

df^ 


< -(- K{Aqn+i - 2Ag„ -I- Ag„_i), 


> -U''{ri)Aqn + K(Ag„+i - 2Ag„ -f Ag„_i) 


respectively and due to reversion symmetry = (y_„ the further analysis can be restricted to n > 0. 
Substituting qn{t) = (—l)"())„/(t) into Eqs. (52) and (53) yields 


< — UJn -\- K 


and 


ll 

df^ 


^>-|C/"(,) + K 


+ 4'n+2 f’n-l + 4>n 

4>n + f’n+l (fn + (fn+l 

4>n-i-l + 4>n+2 , + <t>n 


-h2 


Pn + y'ra+1 (Pn + (fn+l 

with C = min„gz+ {+ 2 } and D = max„gz+ { 


) / < - + kC) f, 

f > - (U'fr,) + kD) f, 

4>n-l-\-4>n _|_ qI 


(52) 


(53) 


(54) 


( 55 ) 






















Thus, by the comparison principle the solutions f{t) are sandwiched as h(t) < f{t) < H{t) where the upper and 
lower bound are given by 


and 


H{t) - i ^ sin(Ubl), 

Ub 

(56) 

H{t) = /(O) cos(Ubi)) 

(57) 

h{t) = i ^ sin(wbt), 

(58) 

Ht) = /(0)cos((Dbt), 

(59) 

and h{t) > 0 for 0 < < < 7r/d)b, where 


^b = U”{r]) -1- kD, 

(60) 


(61) 


The remainder of the proof regarding the time periodicity of the resulting localised solutions proceeds in an analogous 
way as above for Theorem 1. 

Conclusively, spatially localised and time-periodic solutions for out-of-phase motion in hard on-site potentials result 
which satisfy 


\qnit)\ > \qn-i{t)\, n > 0, 
\qn+i{t)\ > \qn{t)\, n < 0 


(62) 

(63) 


for t € K and {qn(t + Th),qn{t + ?&)) = iqn{t), qn{t)) for n S Z with period Tf, = The equality sign in (62) 


and (63) is valid only at the instants of time, tfcgz, when the oscillators pass simultaneously through zero coordinate, 
corresponding to the minimum position of the on-site potential, i.e. q„(tfe) = 0. Using similar arguments as in the 
part of the proof of Theorem 1 that concerns the frequencies of oscillations, Wb, one infers from (61) that breather 
frequencies lie above the upper edge of the linear band fulfilling the relation 


Wfe > 


- 4k. 


(64) 


which concludes the proof. 

□ 


IV. NUMERICAL ANALYSIS 

In this section, we verify, by means of a numerical analysis, the validity of assumptions ( [T^ and ( [4^ . To this 
aim, we calculate breathers by using the numerical methods developed in e.g. m- For examples of soft and hard 
on-site potentials we consider the sine-Gordon potential: U{qn) = 1 —cos( 9 „) and potential: U{qn) = g^/2-)-(;()/4, 
respectively. 

We quantitatively describe the validity of the above assumptions by monitoring the two following coefficients related 
to the relative norms of the exact breather profile ({ 9 n( 0 )}) and its approximation ({<(>»}): 


IVoo = 1- 




ll^nllc 


W = 1- 




n||2 


(65) 


Consequently, the lower the values of N^o and N 2 , the better the approximation made in (15) and (45). 

Figure shows the dependence of the validity coefficients with respect to the frequency Wb and coupling constant 
K for the sine-Gordon potential. Notice that the maximum value of k is attained when 3ujb coincides with the upper 
edge of the phonon band. From this figure, we can observe that, for a given value of the frequency, the minima 
coefficients occurs at the anti-continuum limit (that is k = 0); in addition, the coefficients exhibit a relative maximum 
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at a given value of the coupling that decreases with the frequency. Finally, for a given value of the coupling constant, 
the coefficients decrease when the frequency approaches the phonon band. 

A comparison of the breather profiles with the corresponding approximation is shown in Fig. where two examples 
of two extreme cases (i.e. one for which the coefficients are maximal and another one for which the coefficients are 
minimal) are depicted. There is excellent agreement except at n = 0. 




FIG. 1: The top panels show the dependence of the coefficients Noo and N 2 with respect to the frequency ujb and the coupling 
constant k for the sine-Gordon potential; notice that ujq = 1. An alternative representation of the dependence of the coefficients 
with the coupling constant for four selected frequencies is depicted in the bottom panels. 

For the hard cj)'^ potential, we have a some differences with respect to the sine-Gordon one, as demonstrated in Fig. 

First of all, the validity coefficients possess, for a given value of the frequency, a monotonically decreasing behaviour 
so that the maximum value occurs at the anti-continuum limit and becomes zero (i.e. perfect matching) when the 
breather frequency reaches the phonon band. Secondly, both and N 2 coefficients exhibit the same behaviour. 
Figure]^ shows the profile of two breathers in a similar fashion to Fig. 


V. SUMMARY 

We have proven the existence of time-periodic spatially localised solutions in general infinite nonlinear KG lattice 
systems utilising the comparison principle for differential equations. The time-dependence of the localised solutions is 
of uniform, i.e. site-independent, functional form being valid for solutions near the end points of the spectral phonon 
band. That is, the frequencies ujb of the localised solutions in inhnite KG lattices with soft (hard) on-site potentials 
are close to the lowest (highest) normal mode frequency. To facilitate the existence proof, it is shown that there 
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(a) (b) 



FIG. 2: Breather profiles for the sine-Gordon potential in two extreme cases. Panel (a) corresponds to a breather with ujh = 0.75 
and K = 0.47 whose validity coefficients are N^o = 0.0436 and N 2 = 0.0344. Panel (b) shows a breather with uJb = 0.99 and 
K = 1.95, characterized by N^o = 0.0017 and N 2 = 0.0012. Blue line and dots correspond to the exact profile {q„(0)} and the 
approximated solution {4>n} is depicted by red circles. 


exist spatially localised lattice states whose amplitudes of periodic oscillations at each lattice site are bounded from 
above and below by the oscillating solutions of two auxiliary linear equations. For soft (hard) on-site potentials the 
in-phase (out-of-phase) breathers oscillate with frequencies slightly below (above) the lower (upper) value of the linear 
spectrum of phonon frequencies. With a numerical analysis of the features of solutions with frequencies close to the 
edges of the the linear phonon spectrum we have verified that the ensuing localised solutions are of (virtually) uniform 
time-dependent functional form. 
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FIG. 4: Breather profiles for the sine-Gordon potential in two extreme cases. Panel (a) corresponds to a breather with cjf, = 3 
and K = 0.5 whose validity coefficients are Nao = 0.0340 and N 2 = 0.0339. Panel (b) shows a breather with ujf, = 2 and 
K = 0.7, characterized by Nao = 0.0043 and N 2 = 0.0031. Blue line and dots correspond to the exact profile {9n(0)} and the 
approximated solution {4>n} is depicted by red circles. 





























